Slow relaxation to equipartition in spring-chain systems 
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In this study, one-dimensional systems of masses connected by springs, i.e., spring-chain systems, 
are investigated numerically. The average kinetic energy of chain-end particles of these systems is 
larger than that of other particles, which is similar to the behavior observed for systems made of 
masses connected by rigid links. The energetic motion of the end particles is, however, transient, 
and the system relaxes to thermal equilibrium after a while, where the average kinetic energy of 
each particle is the same, that is, equipartition of energy is achieved. This is in contrast to the case 
of systems made of masses connected by rigid links, where the energetic motion of the end particles 
is observed in equilibrium. The timescale of relaxation estimated by simulation increases rapidly 
with increasing spring constant. The timescale is also estimated using the Boltzmann- Jeans theory 
and is found to be in quite good agreement with that obtained by the simulation. 
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I. INTRODUCTION 

Constrained systems are those having constraints on 
their degrees of freedom. When a constraint is imposed 
on spatial coordinates, the constraint is called a "holo- 
nomic constraint" [J . Constrained systems are useful 
and simple, because of which they are widely used as 
model systems. An example of such a model system is a 
freely jointed chain , which is a model composed of 
one-dimensional chain (ID) of N masses such that the 
distances between adjacent masses are constant. The 
freely jointed chain is known as a simplified model of 
polymers. In computational software packages such as 
CHARMM (I and AMBER 6] for molecular dynamics 
calculations, the SHAKE and RATTLE algorithms en- 
able one to treat model systems as constrained ones by 
setting distances between atoms as constant. As a re- 
sult, the computational task can be made much easier, 
and then, physically important results can be obtained. 
In molecular dynamics calculations, water molecules are 
often treated as having a fixed shape, where the length of 
bonds between hydrogen and oxygen atoms is fixed 
In both cases, we replace the bond between atoms with 
a rigid link when the frequency of bond vibration is ex- 
tremely high. 

It is known that in a constrained system, the equipar- 
tition of energy occurs in a somewhat complicated way, 
and the average kinetic energies of particles {l/2mivf) 
can take different values for particles located at differ- 
ent places in the system P. In the generalized form 
of equipartition [l3| , what is equal among degrees of 
freedom is not {l/2mivf) but (l/2pi|^), where is the 
momentum conjugate to the generalized coordinate qi of 
the I'th degree of freedom, and K is the kinetic energy 
of the system. In constrained systems, K depends on 
coordinates and pi is no longer equal to rmvi; therefore 
l/2pi^, whose average takes the same value for all i, is 
not equal to l/2mivf. 



Recently, we found that for a chain-type system, 
termed a planar chain model, the average kinetic energy 
of each particle differs systematically; that is, particles 
near both the ends of the chain have relatively large av- 
erage kinetic energies Q . This model consists of masses 
connected by rigid links; since the distances between ad- 
jacent masses are fixed in this model, it is a constrained 
system. It is the constraint that induces the nonunifor- 
mity of average kinetic energies. The abovementioned 
energetic motion of end particles observed in the model 
would be useful in understanding the dynamical behavior 
of chain- type systems such as polymers DNAs, pro- 
teins [l^l , and some artificial objects such as manipulator 
arms of spacecraft. 

Thus far, we have described the behavior of con- 
strained systems. However, strictly speaking, a con- 
strained system or rigid link does not exist in the real 
world. When the potentials of a system are somewhat 
steep, i.e., spring constants are quite large and the fre- 
quency of bond vibration is reasonably high, we approxi- 
mate the bonds with rigid links. A rigid link or holonomic 
constraint is an idealized limit of a stiff spring. 

However, if we replace the rigid links with springs, say 
spring-chain model, the usual expression of equipartition 
of energy, i.e., {^rrnvf) = ^ksT, holds regardless of the 
magnitude of the spring constant fc, where D denotes the 
spatial dimension. In other words, although a spring- 
chain system, made by replacing rigid links in the planar 
chain system with springs, appears to behave like the pla- 
nar chain model when the spring constants are large, the 
behavior of energetic motion of end particles cannot be 
reproduced by applying equilibrium statistical mechanics 
to the model. 

Then it would be interesting to know whether the large 
average kinetic energies of end particles observed in the 
planar chain model can also be found in the real world or 
whether it is an artifact observed only in mathematical 
models and is never observed in the real world. If the 



former is true, then we can expect to observe an inter- 
esting feature that the energy distribution of many-body 
systems shows nonuniform behavior. Since it is natural 
to consider that the dynamical behavior of a stiff spring 
is similar to that of a rigid link at least for a finite time, 
the solution to the problem will be the knowledge of the 
relaxation properties to equilibrium. 

With this background, the aim of this study is to ex- 
amine the property of relaxation to equipartition for the 
spring-chain system, particularly for a large value of the 
spring constant k. We measure the relaxation time t^eiax 
and investigate its relation with the spring constant k. 
Further, we estimate t^eiax by the Boltzmann- Jeans the- 
ory and compared it with the value measured. 

This paper is organized as follows. In Sec. II we in- 
troduce the model, the spring-chain system. In Sec. Ill 
we briefly describe the method of numerical computa- 
tion. The results are shown in Sec. IV. The final section 
is devoted to the summary and discussion. 
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II. MODEL 



FIG. 1: A schematic illustration of spring-chain model. 



We now introduce the spring-chain system. It is com- 
posed of N particles (masses) connected by — 1 mass- 
less springs. The masses can rotate smoothly on the xy- 
plane, as shown in Fig. [TJ The system is defined by the 
following Lagrangian L: 
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where rrii is the mass of the «'th particle, = {xi,yi) 
represents the position of the i'th particle, ki and £i are 
the spring constant and natural length of the i'ih spring, 
respectively. U is an external potential. In this paper we 
consider the case oi rrii = m, ki = k, and £i = £ for all 
i. We set m = 1 and £ = 1. The spring-chain model is a 
kind of beads-type models, which are used as models of 
polymer and protein |13l - [l7j . 



III. METHOD OF NUMERICAL INTEGRATION 

We integrate the equation of motion of the model 
by a fourth-order symplectic integrator that is the 
composition of three successive second-order symplec- 
tic integrators. External potential ^] Uil^) = 

is applied in order to 



0.01 Ef=T 



break the rotational symmetry and thus prevent the 
conservation of angular momentum. Here N^aii = 4, 
a = 4:N£, Rj EE (i?,0), (0,i?), (-i?,0), (0,-i?), R = 
N£ + V7? 



N^£^. 

Throughout this paper, the following parameters and 
initial condition are set rrii = m = l, £i = £= l for all 



i and Xi = i — {N + l)/2, yi = for all i, respectively. 
The values of system size and initial momentum for each 
simulation set will be defined in the subsequent sections. 



IV. RESULT 

Next, we briefiy summarize the relation between our 
numerical simulation and thermal equilibrium. Under 
most of the initial conditions considered for the present 
simulations, the system undergoes chaotic motion. Since 
the model has no conserved quantities other than the 
total energy, one may think that the states of the sys- 
tem attained in the course of a long-duration are well 
approximated by a microcanonical distribution. In that 
case, the distribution of the state of each particle in the 
chain can be approximated by a canonical distribution at 
a certain temperature, by considering the other particles 
in the chain as a heat bath. That is, the long-term aver- 
age of kinetic energy of each particle, Ki{t), is equal to 
thermal average (Ki). Then, according to the principle 
of equipartition of energy, the average kinetic energy of 
each particle is the same: 

K-{t)^^l^'^{x,{tr+utr)dt' 

^ {K,) = i / K,exp{-PH)dpdq 

t-s-oo Z J 

^keT, (2) 

where Z = J exp{—(3H)dpdq, (3 = l/fc^T, fc^ is the 
Boltzmann constant, and T is the temperature. Since 
our aim is to investigate the property of relaxation to 
equipartition, we define the following quantity in order 
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FIG. 2; Ki{t) vs. t for fc = 10''. A'' = 8. Plots are measured 
at 10 successive times: t = 10^ (blue circles), 2 x 10^ (red 
squares), 4 x 10^ (yellow diamonds), 8 x 10^ (green triangles), 
and 1.6 x 10® (blue inverted triangles). The time step of 



integration St = 4 x 10 ^. 
i - (AT + l)/2, = 0, pf 
(i>2). 



The initial conditions are Xi — 
= 0, and v\ = 0.1, p\ = -0.1 
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FIG. 3: Time evolution of A(f) for fc = 1 (blue circles), 10 (red 
squares), 10* (yellow diamonds), and 10^ (green triangles). 
The other parameters and initial conditions are the same as 
those mentioned in the caption of Fig. [2] The dashed line 
represents 0.001/t. 



to measure how close the system is to equipartition: 

n 2 



A(t) 



1 ^ 




(3) 



If A(t) = 0, then Ki{t) = Kq for all i. Similar quantities 
have been used to measure the de gree of equipartition 
in studies on a supercooled liquid flSf . self-gravitating 
systems p^| - [2]| . and proteins 

Figure [2] shows the time evolution of the profile of 
{ifi(t), i — \^2, ■ ■ ■ , N — %} . It is clearly observed that 
in the initial stage of the time evolution, the average ki- 
netic energy of all particles is not equal; rather, particles 
near both the ends of the chain have a larger average Ki. 
The profile is similar to that of the rigid link, i.e., the 
planar chain model [1]. Then, as time progresses, differ- 
ences in average Ki among particles gradually decrease 
and tend to zero, and equipartition is achieved. Figure [3] 
shows the relaxation of A(t) [Eq. ([3])] for the data con- 
sidered in Fig. [2j We observe that the system relaxes to 
equilibrium with the progress of time. 

The physical process of relaxation can be understood 
by examining the kinetic energy in greater detail. We 
rewrite the Hamiltonian as 
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FIG. 4: Time evolution of A'^'(t;T) (blue circles) and 
Kyit/Krot (red squares). N ^ 8, K = 100, and St = 10"^. 
r = 2 X 10^. The initial conditions are Xi ^ i — {N + l)/2, 
= 0, pf = 0, = 0.1 {i = 1), and pf = -0.1 (2 < i < 8) . 
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If equipartition is achieved, (Kyu,) = (Krot), because 
the model has the same number of springs and angles. 
Figure 0] shows the temporal evolution of A^^^ {t; r) and 
Kmb{t\T) / Krot{t\T) for k = 200. Here, the time average 
with two arguments, /(t; r), is defined as 
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FIG. 5: k dependence of average relaxation time treiax- At 
each k value, 15 samples are taken. Threshold value Ao — 
10~^. The time step St of integration is scaled as St = 0.125 * 
VW/Vk. The error bars represent standard deviations of 
^r'eiax'''''' for 15 Samples. The dashed line shows treiax ~ 5.52 x 
10* exp(0.415v^) obtained by fitting. 



FIG. 6: Distribution of treiax for k = 100 (blue circles) and 
k = 200 (red squares). For each value of k, the numbers of 
samples are 200. Both distributions are shown on a semi-log 
scale. The dashed lines show P{treiax) = cexp(— atreiai:). 
For k = 100, c = 0.544 and a = 5.39 x 10"^ For k = 200, 
c = 0.271 and a = 1.54 x 10"^ 
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We find that the system relaxed to equipartition on a 
timescale similar to that on which the rotational energy 
Kj-ot transformed into vibrational energy K^ny. 

As mentioned earlier, the aim of this study is to exam- 
ine the property of relaxation to equipartition of energy 
for the model expressed in Eq. ([T]). First, for each sample 
orbit starting from different initial condition, vi^e define 



the relaxation time 



{sample) 
relax 



as the time required for A(t) 



to decay below a critical value Aq . We define the average 
relaxation time over N^ampie orbits as 
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This time is a measure of the relaxation time to equipar- 
tition. 

Figure [5] shows the plot of the dependence of treiax on 
the spring constant k. We observe that as the stiffness of 
the spring increases, the relaxation time increases rapidly. 
That is, systems with hard springs or a steep potential 
show rigid-like behavior of energetic particles near the 
chainends for a very long time, as shown in Fig. [2l before 
relaxing to equipartition. 

Here, we mention a technical detail about the numer- 
ical integration used for obtaining the plot in Fig. [S] 
With increasing spring constant fc, the period of bond- 
stretching vibration decreases in proportion to 1/ Vfc. 
Therefore, for large values of fc, the magnitude of the 
time step of numerical integration should be reduced. We 



confirmed that treiax converges at 5t = 0.125 for k = 10. 
Thus, we used 5t = 0.125 * TTO/v^ for each k. 

On changing the initial conditions, the relaxation 
changes, and we obtain a distribution of 



time t 



[sample) 



{sample) 



relax 



). Figure [6] shows the dis- 
100. For both 
cases, the histograms show exponential decay expressed 
as 



'-relax ' dcuotcd aS Pitrelax 

tribution of treiax for fc = 20 and k 



Pitrelax 

) CX eyiY>{-atrelax) , 



(9) 



which suggests the existence of a characteristic timescale 
for the relaxation. 

Next, we analyze the results of the abovementioned cal- 
culation using the concept of the Boltzmann- Jeans the- 
ory (also known as Boltzmann- Jeans conjecture) [23l - [28l |. 
The essence of this theory is roughly described as follows. 
(For a detailed description of the theory please refer to 
|23l ) . ) Suppose we have a system described by a Hamilto- 
nian, which has two subsystems Hf and H^, their typical 
time scale being r/ and r^., respectively. Here, subscripts 
/ and s denote "fast" and "slow," respectively. Let us 
call Hf and Hs as a "fast subsystem" and "slow sub- 
system," respectively. If the timescales of the fast and 
slow subsystems differ greatly, i.e., Tg/rf ^ 1, then the 
timescale for the occurrence of energy exchange between 
these two subsystems is on the order of 



texch Z exp(c— ) . 



(10) 



That is, energy exchange occurs after a long time. In the 
case of the spring-chain system [Eq. ([1])] with a large 
spring constant k, the fast and slow subsystems corre- 
spond to bond vibration and relative rotation, respec- 
tively. Since the typical timescale of bond vibration is on 
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the order of 2'KyJm/k and that of rotation is assumed to 
be constant, we have 



texch exp(c • 



exp 



(11) 



Since the relaxation to equipartition occurs by energy 
transfer from rotation to vibration (as we observed be- 
fore), we can consider that ie^c/i mentioned above is es- 
sentially the same as the relaxation time treiax- 



relax 



exp 



(c'Vfc) , 



(12) 



Now, we examine whether treiax obtained by the sim- 
ulation obeys Eq. p^ . The result is already shown in 
Fig. [5l hog{treiax) IS proportional to ^/k; therefore the 
interpretation by the Boltzmann- Jeans theory is appro- 
priate. 

This theory can also be used for interpreting the his- 
togram of t^reiax^'^'^ ■ Since treiax IS defined from the aver- 
age of a number of samples, we have 

^relax^i^relax^^^relax ' 



(13) 



where P{t) is the distribution of 4e?a"^''^'' ■ l'^ adopt the 
exponential form for the distribution P [Eq. ([9])], then 



relax 
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(14) 



where a is the coefficient that appears in Eq. 

Combining Eqs. (IT^ and p^ . the relation between 
the coefficient a and the spring constant should be 



a oc exp 



i-c'Vk) . 



(15) 



A comparison between the estimation [Eq. (|15p ] and 
data is shown in Figure [71 from which we find that they 
are in good agreement. Thus, the fact that relaxation to 
equilibrium takes quite a long time to occur can be inter- 
preted as the outcome of the Boltzmann- Jeans theory. 



V. SUMMARY AND DISCUSSION 

In this paper, we have numerically shown the occur- 
rence of energetic motion of end particles for a ID chain 
of point masses connected by hard springs. The timescale 
at which the energetic motion is observed depends on 
the spring constant, and this timescale lengthens with 
increasing the spring constant. Relaxation to equilib- 
rium occurs as rotational kinetic energy is converted into 
vibrational kinetic energy. 

The timescale of relaxation is estimated using the 
Boltzmann-Jeans theory, which describes the energy ex- 
change rate in a system in which fast and slow motions 
coexist. In the case of our model, fast motion corresponds 
to bond vibration by a hard spring and slow motion cor- 
responds to rotation and deformation of the chain. The 
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FIG. 7: Plot of a vs k, obtained from 300 samples. The 
dashed line shows logo — ci — C2y/k, where ci = —3.93 and 
C2 = 0.124. 



result of our numerics is that the timescale is estimated as 
an exponential of the square root of the spring constant, 
which coincides well with simulation data. 

The energetic motion of end particles is also observed 
in other systems , e.g., a planar chain system Q and mul- 
tiple pendulum |29l - [3ll | . These systems consist of masses 
connected by rigid links, and therefore, they are con- 
strained systems. If we consider a rigid link as a lim- 
iting case of a hard spring when the spring constant k 
is large (fc — > oo), then the spring-chain system exam- 
ined in this study becomes a planar chain system when 
k — ^ oo. Therefore, it is natural that the spatial energy 
distribution of the spring-chain system considered in this 
study resembles those of systems with rigid links. 

However, the equilibrium behaviors of spring-chain 
system and rigid- link chain system are quite different. 
In the spring-chain system, equipartition achieved is of 
the usual form; that is, the average kinetic energy of 
each particle is the same, regardless of the magnitude 
of k. Then, chain-end particles in thermal equilibrium 
should not exhibit the energetic behavior. In contrast, in 
the rigid-link-chain system, chain-end particles behave 
energetically even in thermal equilibrium Q. The ener- 
getic behavior of end particles in the spring-chain sys- 
tem is a transient behavior before the system relaxes to 
thermal equilibrium. What is important here is that the 
relaxation time increases with increasing the spring con- 
stant k and that the relaxation time eventually diverges 
at fc ^ oo. Therefore, a large k value provides a good 
opportunity to observe energetic behavior in spring-chain 
systems. 

The relaxation time to equipartition changes if we 
change the initial conditions. We found that the the 
relaxation time is distributed according to the exponen- 
tial form for large relaxation time. This implies that 
relaxation occurs almost randomly. That is, the sys- 
tem moves on the energy surface in a random way and 
happens to encounter at which the system can divert to- 
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ward the equipartition state. This situation is in contrast 
to the process of slow dynamics often observed in many 
Hamiltonian systems, where the relaxation time is often 
distributed according to the power law. 

Slow relaxation is often observed in Hamiltonian sys- 
tems, and in most cases, it is accompanied by strong 
temporal correlation caused by sticky or stagnant mo- 
tion around KAM tori (regular orbits) and their rem- 
nants and l//-type fluctuations. Such slow relaxation 
is observed in area-preserving mappings and some other 
high-dimensional systems [32 40]. Such slow relaxation 
is explained by a hierarchical structure generally found 
in nearly integrable Hamiltonian systems. In this sense, 
it is quite common to find slow relaxation in Hamiltonian 
systems, if the Hamiltonian is nearly integrable. 

However, the spring-chain model considered in this 
study is not nearly integrable. In addition, the systems 
which show slow relaxation in real world are not always 
nearly integrable. In contrast, the Boltzmann- Jeans the- 
ory is still applicable to Hamiltonian systems that are 
not nearly integrable and is able to explain the slow re- 
laxation. Similar to the research by Shudo et al. [4l[, 
the result of our study also implies that the Boltzmann- 
Jeans theory can describe slow relaxation in more general 
systems other than the nearly integrable ones. 



We used the Boltzmann- Jeans theory [23l - l28| to es- 
timate the relaxation time. Although the concept of 
the theory dates back to the 19th century [43 - |44| . its 
importance is not very familiar and not many exam- 
ples of the application of this theory have been demon- 
strated [4ll . l45l - l49| . Because the results of this study show 
that there is good agreement between numerical data and 
theoretical estimation, this study can be a good example 
of the applicability of the Boltzmann- Jeans theory. 

In this study, we showed that chain-end particles be- 
have energetically even for systems with a finite but large 
spring constant. This implies that similar behavior can 
be observed for natural chain-type systems that are not 
made of rigid links but whose intrachain potential be- 
tween elements is very steep. The results of this study 
are expected to have many useful applications to poly- 
mers, proteins and some artificial objects. 
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